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Abstract 



Vacuum expectation values of the field square and the energy-momentum tensor for the 
electromagnetic field are investigated for the geometry of a wedge with a coaxal cylindrical 
boundary. All boundaries are assumed to be perfectly conducting and both regions inside 
and outside the shell are considered. By using the generalized Abel-Plana formula, the 
vacuum expectation values are presented in the form of the sum of two terms. The first one 
corresponds to the geometry of the wedge without the cylindrical shell and the second term 
is induced by the presence of the shell. The vacuum energy density induced by the shell 
■ is negative for the interior region and is positive for the exterior region. The asymptotic 

behavior of the vacuum expectation values are investigated in various limiting cases. It is 
shown that the vacuum forces acting on the wedge sides due to the presence of the cylindrical 
boundary are always attractive. 



PACS numbers: 03.70.+k 

O 

! 1 Introduction 

The Casimir effect is among the most interesting macroscopic manifestations of vacuum fluctu- 
ations. It have important implications on all scales, from cosmological to subnuclear, and has 
become in recent decades an increasingly popular topic in quantum field theory. In addition to 
its fundamental interest the Casimir effect also plays an important role in the fabrication and 
operation of nano- and micro-scale mechanical systems. The imposition of boundary conditions 
on a quantum field leads to the modification of the spectrum for the zero-point fluctuations and 
results in the shift in the vacuum expectation values for physical quantities such as the energy 
density and stresses. In particular, the confinement of quantum fluctuations causes forces that 
act on constraining boundaries. The particular features of the resulting vacuum forces depend 
on the nature of the quantum field, the type of spacetime manifold, the boundary geometries 
and the specific boundary conditions imposed on the field. Since the original work by Casimir [1] 
many theoretical and experimental works have been done on this problem (see, e.g., [2j [31 HI [5] 
and references therein). Many different approaches have been used: mode summation method 
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with combination of the zeta function regularization technique, Green function formalism, mul- 
tiple scattering expansions, heat-kernel series, etc. Advanced field-theoretical methods have 
been developed for Casimir calculations during the past years [Sid [5]. However, there are still 
difficulties in both interpretation and renormalization of the Casimir effect. Straightforward 
computations of geometry dependencies are conceptually complicated, since relevant informa- 
tion is subtly encoded in the fluctuations spectrum [8\. Analytic solutions can usually be found 
only for highly symmetric geometries including planar, spherically and cylindrically symmetric 
boundaries. Recently the Casimir energy has been evaluated exactly for several less symmetric 
configurations of experimental interest. These include a sphere in front of a plane and a cylinder 
in front of a plane [9]. 

A great deal of attention received the investigations of quantum effects for cylindrical bound- 
aries. In addition to traditional problems of quantum electrodynamics under the presence of 
material boundaries, the Casimir effect for cylindrical geometries can also be important to the 
flux tube models of confinement |10[ [TT] and for determining the structure of the vacuum state 
in interacting field theories [12] . The calculation of the vacuum energy of electromagnetic field 
with boundary conditions defined on a cylinder turned out to be technically a more involved 
problem than the analogous one for a sphere. First the Casimir energy of an infinite perfectly 
conducting cylindrical shell has been calculated in Ref. [13J by introducing ultraviolet cutoff 
and later the corresponding result was derived by zeta function technique [14L I15|, [T6] . The local 
characteristics of the corresponding electromagnetic vacuum such as energy density and vacuum 
stresses are considered in [17] for the interior and exterior regions of a conducting cylindrical 
shell, and in |18| for the region between two coaxial shells (see also [19J). The vacuum forces 
acting on the boundaries in the geometry of two cylinders are also considered in Refs. |20| . 
The scalar Casimir densities for a single and two coaxial cylindrical shells with Robin boundary 
conditions are investigated in Refs. [214 122]. Less symmetric configuration of two eccentric per- 
fectly conducting cylinders is considered in Ref. [23]. Vacuum energy for a perfectly conducting 
cylinder of elliptical section is evaluated in Ref. [24] by the mode summation method, using the 
ellipticity as a perturbation parameter. The Casimir forces acting on two parallel plates inside 
a conducting cylindrical shell are investigated in Ref. |25| . 

Aside from their own theoretical and experimental interest, the exactly solvable problems 
with this type of boundaries are useful for testing the validity of various approximations used 
to deal with more complicated geometries. From this point of view the wedge with a coaxial 
cylindrical boundary is an interesting system, since the geometry is nontrivial and it includes 
two dynamical parameters, radius of the cylindrical shell and opening angle of the wedge. This 
geometry is also interesting from the point of view of general analysis for surface divergences 
in the expectation values of local physical observables for boundaries with discontinuities. The 
nonsmoothness of the boundary generates additional contributions to the heat kernel coefficients 
(see, for instance, the discussion in |26[ l27l [28] and references therein). The present paper is 
concerned with local analysis of the vacuum of the electromagnetic field constrained to satisfy 
perfectly conducting boundary conditions on boundary surfaces of a wedge with a coaxial cylin- 
drical boundary. Namely, we will study the vacuum expectation values of the field squares and 
the energy-momentum tensor for the electromagnetic field for both regions inside and outside 
the cylindrical shell. In addition to describing the physical structure of the quantum field at 
a given point, the energy- momentum tensor acts as the source of gravity in the Einstein equa- 
tions. It therefore plays an important role in modelling a self-consistent dynamics involving the 
gravitational field. The vacuum expectation value of the square of the electric field determines 
the electromagnetic force on a neutral polarizable particle. Some most relevant investigations to 
the present paper are contained in Refs. [2] [291 USB Ell E21 [33] , where the geometry of a wedge 
without a cylindrical boundary is considered for a conformally coupled scalar and electromag- 
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netic fields in a four dimensional spacetime. The total Casimir energy of a semi-circular infinite 
cylindrical shell with perfectly conducting walls is considered in [M] by using the zeta function 
technique. For a scalar field with an arbitrary curvature coupling parameter the Wightman 
function, the vacuum expectation values of the field square and the energy-momentum tensor 
in the geometry of a wedge with an arbitrary opening angle and with a cylindrical boundary 
are evaluated in |35|, 136] . Note that, unlike the case of conformally coupled fields, for a general 
coupling the vacuum energy-momentum tensor is angle-dependent and diverges on the wedge 
sides. Our method here employs the mode summation and is based on a variant of the gen- 
eralized Abel-Plana formula [37] (see also Refs. |19l I38j ). This enables us to extract from the 
vacuum expectation values the parts due to a wedge without the cylindrical shell and to present 
the parts induced by the shell in terms of strongly convergent integrals. Note that the closely 
related problem of the vacuum densities induced by a cylindrical boundary in the geometry of 
a cosmic string is investigated in Refs. [39, 3UJ for both scalar and electromagnetic fields. 

We have organized the paper as follows. In the next section we describe the structure of 
the modes for a wedge with a cylindrical shell in the region inside the shell. By applying 
to the corresponding mode sums the generalized Abel-Plana formula, we evaluate the vacuum 
expectation values of the electric and magnetic field square. Various limiting cases of the general 
formulae are discussed. Section [3] is devoted to the investigation of the vacuum expectation 
values for the energy-momentum tensor of the electromagnetic field in the region inside the 
shell. The additional vacuum forces acting on the wedge sides due to the presence of the 
cylindrical boundary are evaluated. In section [5] we consider the vacuum densities for a wedge 
with the cylindrical shell in the exterior region with respect to the shell. Formulae for the 
shell contributions are derived and the corresponding surface divergences are investigated. The 
vacuum forces acting on the wedge sides are discussed. The main results are summarized and 
discussed in section [5] 

2 Vacuum expectation values of the field square inside a cylin- 
drical shell 

Consider a wedge with the opening angle (j)Q and with a coaxial cylindrical boundary of radius 
a (see figure CO) assuming that all boundaries are perfectly conducting. In accordance with the 
problem symmetry, in the discussion below the cylindrical coordinates (r, <j), z) will be used. We 
are interested in the vacuum expectation values (VEVs) of the field square and the energy- 
momentum tensor for the electromagnetic field. Expanding the field operator in terms of the 
creation and annihilation operators and using the commutation relations, the VEV for a quantity 
F {Ai, A^} bilinear in the field can be presented in the form of the mode-sum 



where {A a i,A* ak } is a complete set of solutions of the classical field equations satisfying the 
boundary conditions on the bounding surfaces and specified by a set of quantum numbers a. 

In accordance with formula (pQ), for the evaluation of the VEVs for the square of the electric 
and magnetic fields and the energy-momentum tensor, the corresponding eigenfunctions are 
needed. In this section we consider the region inside the cylindrical shell (region I in figure 
[1]). For the geometry under consideration there are two different types of the eigenfunctions 
corresponding to the transverse magnetic (TM) and transverse electric (TE) waves. In the 
discussion below we will specify these modes by the index A = and A = 1 for the TM and TE 
waves respectively. In the Coulomb gauge, the vector potentials for the TM and TE modes are 




(1) 
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Figure 1: Geometry of a wedge with a coaxial cylindrical boundary with radius a. 

given by the formulae 

A = 3 I t 1 /^) ( 72es + ikVt ) J i\ m \ ^ sin (<? m ^) exp [i (kz - cut)] , A = 

a \ -e 3 xV t {j qlml (jr)cos(qm(j))exp[i(kz-u;t)]}, A = l ' 

where e3 is the unit vector along the axis of the wedge, Vt is the part of the nabla operator 
transverse to this axis, J v {x) is the Bessel function of the first kind, and 

u; 2 = 1 2 + k 2 , q = ir/fo. (3) 

In Eq. ([2|), m = 1, 2, . . . for A = and m = 0, 1, 2, . . . for A = 1. The normalization coefficient 
P a is found from the orthonormalization condition for the vector potential: 

2-7T 

dVA a • A;, = — S aa >, (4) 

where the integration goes over the region inside the shell. From this condition, by using the 
standard integral involving the square of the Bessel function, one finds 

2 _ 4qT qm (ja) x j 1/2, m = . . 

irujaj [ 1, m f 

where we have introduced the notation 

T v {x) = x \j' 2 {x) + (1 - u 2 /x 2 )J 2 (x)j ^ . (6) 

Eigenfunctions ([2]) satisfy the standard boundary conditions for the electric and magnetic 
fields, n x E = and n • B = 0, on the wedge sides corresponding to (j) = and 4> = <f)Q, with n 
being the normal to the boundary. The eigenvalues for the quantum number 7 are determined by 
the boundary conditions on the cylindrical shell. From the latter it follows that these eigenvalues 
are solutions of the equation 

4%(ia) = 0, A = 0,1, (7) 

where we use the notations Ju°\x) = J v (x) and Jp (x) = J' u {x). We will denote the corre- 
sponding eigenmodes by 7a = jm)n, n = 1, 2, . . ., assuming that the zeros jm]n are arranged in 
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ascending order. Consequently, the eigenfunctions are specified by the set of quantum numbers 
a = (k, m, A, n). 

First we consider the VEVs of the squares of the electric and magnetic fields inside the shell. 
Substituting the eigenfunctions ([2]) into the corresponding mode-sum formula, we find 

A a ~ , r +0 ° ^ ^ 7 -W 3 T (j W ) 

\ / V Jm,n ^qrn\Jm,n) 



nu ™_n J-oo x_n 1 1 



m=0 J ~°° A=0,1 n=l yj jgg + fc2 a 2 

x^^gO^J^^^r/a)], (8) 

where F = E,B with = A, t/ba = 1 — A, and the prime in the summation over m means that 
the term m = should be halved. In formula (JE]), for a given function f(x), we have introduced 
the notations 

[$(</>),/(*)] = (fcV/x 2 )[$ 2 (0)/' 2 (x) + $' 2 (0)/ 2 (x)/x 2 ]+$ 2 (0)/ 2 (x), (9) 
<7 (1) [$(0),/(*)] = (l + A; 2 r 2 /x 2 )[<I> 2 (0)/ /2 (x) + <I> /2 (0)/ 2 (x)/x 2 ], (10) 



and 



" W) 1 cos(W>), A = l ' 1 j 



The expressions (|8j) corresponding to the electric and magnetic fields are divergent. They may 
be regularized introducing a cutoff function ip^u)) with the cutting parameter \l which makes 
the divergent expressions finite and satisfies the condition ipu{u) — ► 1 for fj, — ► 0. After the 
renormalization the cutoff function is removed by taking the limit [i — ► 0. An alternative way 
is to consider the product of the fields at different spacetime points and to take the coincidence 
limit after the subtraction of the corresponding Minkowskian part. Our approach here follows 
the first method. 

As we do not know the explicit expressions for the zeros jm,n as functions on m and n, and 
the summand in formula ([8]) is strongly oscillating function for large values of m and n, this 
formula is not convenient for the further evaluation of the VEVs of the field square. In order to 
obtain an alternative representation, we apply to the series over n the generalized Abel-Plana 
summation formula [37] (see also [55] ) 



00 I />oo 

Y. T im{&)f{&n) = 2/ dxf(x) + -Resf( Zj (V) 

n=l ^° Z Jqm(z) 



2vr J t W 

■ 1 qm 



X 



e- gm7vi f(e ni/2 x) + e qmni f{e-™ /2 x) , (12) 



where Y u (z) is the Neumann function and I v (z), K u (z) are the modified Bessel functions. As 
it can be seen, for points away from the shell the contribution to the VEVs coming from the 
second integral term on the right-hand side of (112j) is finite in the limit [i — > and, hence, the 
cutoff function in this term can be safely removed. As a result the VEVs can be written in the 
form 

(0|F 2 |0) = (0 w |F 2 |0 w ) + (F 2 > cyl , (13) 
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where 



\F 2 \0 V 



m=0 



dk 



7 3 ^(w) 



o \/7 2 + k 2 



1 + 



2k 



m 



+ Jqm(nfr) 



-(-1) VF1 co$(2qm<t)) 



' qm 



1 + 



q 2 m 2 

^y2 y»2 



J qm(l r ) 



and 



(^ 2 )cyl 



7T 



2 Z-^ 



A=0,1 



" 3 i^S(xa) *g£(0), I gm (xr)] 



\fx 2 ~ 



k 2 



m=0 " u A=0,1 " ™ -""J™ 1 

In formula Eq. (|15() we have introduced the notations 

G< > [A;, /(x)] = (k 2 r 2 /x 2 ) [<S> 2 (<fr)f 2 (x) + & 2 (0)f 2 (x)/x 2 } + ^> 2 (<fr)f 2 (x) 
= (feV/x 2 -l)[$ 2 (0)/ /2 (x) + $ /2 (0)/ 2 ( a; )/x 2 ]. 



(14) 



(15) 



(16) 
(17) 



oo and the first 
> corresponds to 



The second term on the right-hand side of Eq. (|13p vanishes in the limit a - 
one does not depend on a. Thus, we can conclude that the term (0 W |F 2 |0 V 
the part in the VEVs when the cylindrical shell is absent with the corresponding vacuum state 
|0 W ). Hence, the application of the generalized Abel-Plana formula enables us to extract from 
the VEVs the parts induced by the cylindrical shell without specifying the cutoff function. In 
addition, these parts are presented in terms of the exponentially convergent integrals. 

First, let us concentrate on the part corresponding to the wedge without a cylindrical shell. 
First of all we note that in Eq. (|14p the part which does not depend on the angular coordinate (j) 
is the same as in the corresponding problem of the cosmic string geometry with the angle deficit 
2-7T — (fro (see Ref. [10]), which we will denote by (0 S |F 2 |0 S ). For this part we have 



(0 S |F 2 |0 S 



m=0 



dk 



d-i 



vV + k 2 



1 + 



7^2 



Jqmiir) + 



7 2 m 2 



^2 r 2 "Qm 



J qm (7 r ) 



+ JqmM 



(OmI^IOm) 



(g 2 - l)(g 2 + 11) 

180vrr 4 



(18) 



where (0m|-^ 2 |0m) is the part corresponding to the Minkowskian spacetime without boundaries 
and in the last expression we have removed the cutoff. To evaluate the part in (I14p which 
depends on <fr, we firstly consider the case when the parameter q is an integer. In this case the 
summation over m can be done by using the formula [411 132] 



£' 

m=0 



q-1 



cos(2qm<fr)J 2 (y) 



— J2 Jo(2ysin(0 + 0oO). 
q 1=0 



(19) 



The formulae for the other series entering in Eq. ()14[) are obtained from (|19p taking the deriva- 
tives with respect to <f> and y. In particular, for the combination appearing in the angle-dependent 
part we obtain 



oo 

£' 

m=0 



cos{2qm<fr) 



Jqm(.V) 



1 + 



7 2 m 2 



Jqmiv) 



q-1 

J^jJ(2ysin(0 + ^oZ)). 

1=0 



(20) 
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Substituting this in formula (|14|) . the integrals remained are evaluated by introducing polar 
coordinates in the (k, 7)-plane. In this way one finds 



3(-l) 



VFl 



q-1 



(0 w |i^|0 w ) = (0 S |F 2 |0 S ) ~ " v 4 ^ 4 ^sin- 4 (0 + /7r/g). (21) 

1=0 

The sum on the right hand-side of this formula is evaluated by the double differentiation of the 
relation [H] 

g-l 

^2 cos~ 2 (x + lit/q) = q 2 sin~ 2 (qx + qir/2). (22) 
z=o 

Finally, for the renormalised VEVs of the field square in the geometry of a wedge without a 
cylindrical boundary we find 



(g2_l)(g2 + H) (-1)^ \ 3q 2 , 1 _ 2 



/ w,ren 



1807rr 4 2vrr 4 sin 2 (g0) 



2 sin (q<p) 



(23) 



with t)e\ = 1 and 77^1 = 0. Though we have derived this formula for integer values of the 
parameter q, by the analytic continuation it is valid for non-integer values of this parameter as 
well. The expression on the right of formula f)23|) is invariant under the replacement eft — > (j)Q — <j> 
and, as we could expect, the VEVs are symmetric with respect to the half-plane <f> = 4>o/2- 
Formula f)23|) for F = E was derived in Ref. [32] within the framework of Schwinger's source 
theory. For q = 1 from formula (|23|) as a special case we obtain the renormalized VEVs of the 
field square for a conducting plate. In this case x = rsin0 is the distance from the plate and 
one has 

<**>*-- (24 » 

Another special case q = 1/2 corresponds to the geometry of a half-plane. In (|23p taking the 
limit r — * 00, with xq = rcj)Q being fixed, we obtain the corresponding results in the region 
between two parallel plates located at the points x = and x = xq: 



(F 2, (-l)^TT 3 
\* /2pl,ren ■ 



180xq 2xq sin (ttx/xq) 



2 sin (ttx/xq) 



1 



(25) 



Now, we turn to the investigation of the parts in the VEVs of the field square induced by 
the cylindrical boundary and given by formula (|15p . By using the formula 



x ' dkk m r 'dx x/(x) = ^ r ( m ^) r dxx m+i f(x), (26) 



V* 2 - k 2 2r(f + l 
these parts are presented in the form 

TT 



E E / dxx"^^G^\^)J qm {xr)]. (27) 

m=0 A=0,1 *<ltn \ xa ) 



Here, for given functions f{x) and <&((/)), we have introduced the notations 

G< > /(x)] = $ 2 (0)/' 2 (x) + $ ,2 (0)/ 2 (x)/x 2 + 2$ 2 (0)/ 2 (z), (28) 
G« /(*)] = -^ 2 {<j ) )f 2 (x)-^ 2 (<l ) )f 2 (x)/x 2 . (29) 

As we see the parts in the VEVs induced by the cylindrical shell are symmetric with respect to 
the half-plane <f> = 4>q/2. 
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The expression in the right-hand side of (|27p is finite for < r < a including the points on the 
wedge sides, and diverges on the shell. To find the leading term in the corresponding asymptotic 
expansion, we note that near the shell the main contribution comes from large values of m. By 
using the uniform asymptotic expansions of the modified Bessel functions (see, for instance, 
for large values of the order, up to the leading order, for the points a — r <C a\ sin <j>\ , a\ sin((/> - 
we find 

3(_1)'?fi 



For the points near the edges (r = a, <j) = 0, 0o) the leading terms in the corresponding asymp- 
totic expansions are the same as for the geometry of a wedge with the opening angle (fro = vr/2. 
The leading terms given by formula (|30p are the same as for the geometry of a single plate (see 
(|24p ). They do not depend on </>o and have opposite signs for the electric and magnetic fields. 
In particular, the leading terms are cancelled in the evaluation of the vacuum energy density. 
Surface divergences originate in the unphysical nature of perfect conductor boundary conditions 
and are well-known in quantum field theory with boundaries. In reality the expectation values 
will attain a limiting value on the conductor surface, which will depend on the molecular de- 
tails of the conductor. From the formulae given above it follows that the main contribution to 



(F }cyi are due to the frequencies u < (a 



Hence, we expect that formula (|27p is valid 



for real conductors up to distances r for which (a — r) <C u>q, with ojq being the characteristic 
frequency, such that for to > ojq the conditions for perfect conductivity fail. 

Near the edge r = 0, assuming that r/a <C 1, the asymptotic behavior of the part induced in 
the VEVs of the field square by the cylindrical shell depends on the parameter q. For q > l+rjpi, 
the dominant contribution comes from the lowest mode m = and to the leading order one has 



(^ 2 )c y l 



-1 



\VF1 



2 1 - r iFi 



r \2»?fi 
2a) 



.I h{x) 



(31) 



0, whereas (E 2 



In this case the quantity (B 2 ) cy \ takes a finite limiting value on the edge r 
vanishes as r 2 . For q < 1 + rjpi the main contribution comes form the mode with m 
the shell-induced parts diverge on the edge r = 0. The leading terms are given by the formula 



cyl 

1 and 



<^ 2 > cy i 



(- l)^i g ( r /a) 
2 2 i- l irT 2 (q)a 4 



2(g-l) roo 



dxx 2q+1 



K q (x) K' q {x) 
I q (x) I'Jx) 



(32) 



As for the points near the shell, here the leading divergences in the VEVs of the electric and 
magnetic fields are cancelled in the evaluation of the vacuum energy density. For q = 1 + rjpi the 
main contribution comes from the modes m = 0, 1 and the corresponding asymptotic behavior 
is obtained by summing the right-hand sides of Eqs. (|3ip and (|32p . In accordance with (|23p . 
near the edge r = the total VEV is dominated by the part coming from the wedge without 
the cylindrical shell. Here we have considered the VEVs for the field square. The VEVs for the 
bilinear products of the fields at different spacetime points may be evaluated in a similar way. 

Now, we turn to the investigation of the behavior of the VEVs induced by the cylindrical 
boundary in the limit q S> 1. In this limit the order of the modified Bessel functions is large for 
m ^ 0. By using the corresponding asymptotic formulae it can be seen that the contribution of 
these terms is suppressed by the factor exp[— 2qm ln(a/r)]. As a result, the main contribution 
comes from the lowest mode m = and the VEVs induced by the cylindrical shell are propor- 
tional to q. Note that in this limit the part corresponding to the wedge without the cylindrical 
shell behaves as q 4 . 
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3 Vacuum energy- momentum tensor inside the cylindrical shell 

Now let us consider the VEV of the energy-momentum tensor in the region inside the cylindrical 
shell. Substituting the eigenfunctions ([2]) into the corresponding mode-sum formula, for the non- 
zero components we obtain (no summation over i) 

m\o) = A-sit f +0 ° dk E E (33) 

Z7T a m =0 J ~°° A=0,1 n=l J j + k 2 a 2 



(01^10) = -^-Y! m*n{2qm4) j dk £ (-1 



m=0 " A=0,1 



oo .(A) ^ / .(A) \ 

e J ;;; m(Jm,n) 4(&/"). (34) 



where we have introduced the notations 



/(*)] = (-ir(2A; 2 /7 2 + l) [^(<A)/' 2 (x) + & 2 (<P)f 2 (x)/y 2 } + <5> 2 (<P) f 2 (x) , (35) 
= (-l)^ 2 (^/ /2 (x)-[d> 2 (^ + (-l)^' 2 (0)/x 2 ]/ 2 (x), (36) 

with j = 0, 3 and I = 1,2. As in the case of the field square, in formulae (|33[) and (|34p we 
introduce a cutoff function and apply formula (|12|) for the summation over n. This enables us 
to present the vacuum energy-momentum tensor in the form of the sum 

<0|:^|0) = <0 w |^f |0 W > + (iXb (37) 

where (0 w |jf|0 w ) is the part corresponding to the geometry of a wedge without a cylindrical 
boundary and (T^) cy \ is induced by the cylindrical shell. By taking into account ([26]) . the latter 
may be written in the form (no summation over i) 

= AE' E / dxx^^^F^MJ^xr)], (38) 



m=0 A=0,1 ' 



m=0 A=0,1 170 



dxx Iq m {xr), (39) 

igm [xa) 



with the notations 



/(*/)] = $ 2 (0)/ 2 (y), i = 0,3, (40) 
F» [$(£), /(?/)] = -(-l) i * 2 (0)/' 2 (i/)-[* 2 (0)-(-l) i * ,2 (0)/y 2 ]/ 2 (i/) ) i = l ) 2. (41) 

The diagonal components are symmetric with respect to the half-plane <fi = (f>o/2, whereas the 
off-diagonal component is an odd function under the replacement (j) — > (po — 4>. As it can be 
easily checked, the tensor {Tf) cy \ is traceless and satisfies the covariant continuity equation 
Vfc(T i fc ) cy i = 0. For the geometry under consideration the latter leads to the relations 

|:(^M+r^(2?Wi = 0, (42) 
^ r [r{TlU)+r^-{T 2 ) cyl = (T 2 ) cyl . (43) 
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As it is seen from formula (|39|) . the off-diagonal component (T 2 1 ) cy i vanishes at the wedge sides 
and for these points the VEV of the energy-momentum tensor is diagonal. By using the in- 
equalities I' u {x) < yl + v 2 j x 2 I v (x) and —K' v {x) > a/1 + v 2 j x 2 K v {x) for the modified Bessel 
functions, it can be seen that K' v {x) / I' v {x) + K v (x)/I v (x) < 0. From this relation it follows 
that the vacuum energy density induced by the cylindrical shell in the interior region is always 
negative. 

The renormalized VEV of the energy-momentum tensor for the geometry without the cylin- 
drical shell is obtained by using the corresponding formulae for the field square. For the corre- 
sponding energy density one finds 

(7o)w,ren = g~ ((-E 2 )w,ren + (-B 2 )w,rcn) = -— ^Orr 2 ^ "' ^ 

As we see the parts in the VEVs of the field square which diverge on the wedge sides cancel 
out and the corresponding energy density is finite everywhere except the edge. Formula (I44h 
coincides with the corresponding result for the geometry of the cosmic string (see |44} I45j) with 
the angle deficit 2ir — 4>q and in the corresponding formula q = 2Tr/<j)Q. Other components are 
found from the tracelessness condition and the continuity equation and one has [29j [30] (see also 

(^)w,rcn = - (ga ~Jffi ^ diagfc 1, -3, 1). (45) 

As we could expect this VEV vanishes for the geometry of a single plate corresponding to q = 1. 
In the limit r —* oo, for fixed values xo = rtfo, from (|45p the standard result for the geometry of 
two parallel conducting plates is obtained. 

The force acting on the wedge sides is determined by the component (T| ) rcn of the vacuum 
energy-momentum tensor evaluated for <f> = and 4> = <Po- On the base of formula (I37j) for the 



corresponding effective pressure one has 

P2 = -{Ti)ren\^,=0,<l> a = P2w + P2cyh (46) 

where 

( g 2 -l)(g 2 + ll) 

P2w = 2407V4 ■ ( 47 ) 

is the normal force acting per unit surface of the wedge for the case without a cylindrical 
boundary and the additional term 

p 2cy i = -<T 2 2 ) cyl | 0=o , 0o =-^E' E / dxx * lx), F& 9 m(xr)), (48) 

m=0 A=0,1 ^<i m ( Xa ' 



F w lf{)] _{ ^f 2 (y)/y 2 , a = o 



with the notation 

P(A)rf ^ 1 ~ <1 ' -f' 2 (y)-f(y), 

is induced by the cylindrical shell. From formula (|47|) we see that the corresponding vacuum 
forces are attractive for q > 1 and repulsive for q < 1. In particular, the equilibrium position 
corresponding to the geometry of a single plate (q = 1) is unstable. As regards to the part 
induced by the cylindrical shell, from (|4'8j) it follows that P2 Cy \ < and, hence, the corresponding 
forces are always attractive. 

Now, let us discuss the behavior of the boundary-induced part in the VEV of the energy- 
momentum tensor in the asymptotic region of the parameters. Near the cylindrical shell the 
main contribution comes from large values of m. Thus, using the uniform asymptotic expansions 
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for the modified Bessel functions for large values of the order, up to the leading order, for the 
points a — r <C a\ sin 0| , a\ sin(0o — 4>)\ we find 



( T 0°)cyl 



T^^cyl 



(a 



60vr 2 ( 



( T l)cyl 



(a — r) ' 
60vr 2 a 2 



(50) 



These leading terms are the same as those for a cylindrical shell when the wedge is absent. For 
the points near the edges (r = a, (f> = 0, ^o) the leading terms in the corresponding asymptotic 
expansions are the same as for the geometry of a wedge with the opening angle = vr/2. The 
latter are given by (|45p with q = 2. Near the edge, r — > 0, for the components (no summation 
over i) (T?) cy i, i = 0, 3, the main contribution comes from the mode m = and we find 



cyl 



4vr 2 a 4 



' 3^0*0 
dxx 



-0.0590- 



0,3. 



(51) 



For the components (no summation over i) (T-) CJ \, i = 1,2, when g > 1 the main contribution 
again comes form m = term and one has {T^) cy \ ~ — (Tq)^, i = 1,2. For g < 1 the main 
contribution into the components (T?) cy i, i = 1, 2, comes from the term m = 1 and we have (no 
summation over i) 



i T i)cy\ 



-l)^cos(2#) /r\2(g-i) 



2 2 <2+ivr 2 r 2 (g)a 4 Vo 



dx x 



2q+l 



K' q {x) 
I q {x) I'Jx) 



K q {x) 



1,2. 



(52) 



In this case the radial and azimuthal stresses induced by the cylindrical shell diverge on the 
edge r = 0. In the case q = 1 the sum of the contributions of the terms with m = and m = 1 
given by formulae (|51|) and (|52|) should be taken. For the off-diagonal component the main 
contribution comes from the m = 1 mode with the leading term 



( r 2)cyl 



gsm 



2 2 <2+ivr 2 r 2 (g)a 3 



r \ 2g-l 

a J 



dxx 2q+l 



K q (x) K' q (x) 
Ux) I'Jx) 



(53) 



and this component vanishes on the edge for q > 1/2. 

In the limit g > 1, the contribution of the modes with m ^ 1 is suppressed by the factor 
exp[— 2gmln(a/r)] and the main contribution comes from the m = mode. The leading terms 
are given by the formulae (no summation over i) 



( T i)cy\ 
(TT)cyl 



4vr 2 a 4 .,o 



4vr 2 a 4 



dx x s K T ^ X } I$(xr/a), i = 0, 3, 



(54) 



dxx 3 -^ [/ 2 (W«) + (-l)*/?(xr/a)] , i = 1,2. (55) 



Though in this limit the vacuum densities are large, due to the factor 1/q in the spatial volume, 
the corresponding global quantities tend to finite value. In particular, as it follows from Eq. 
(|55p . in the limit under consideration one has (T/) cy i > 0. Note that in the same limit the parts 
corresponding to the wedge without the cylindrical shell behave as g 4 and, hence, for points not 
too close to the shell these parts dominate in the VEVs. 

In figures [2][5] we have plotted the parts in the VEVs of the energy-momentum tensor induced 
by the cylindrical shell, a 4 (Tf) cy i, as functions of x = (r/a) cos 4> and y = (r/a) sin^>, for a wedge 
with the opening angle </>o = tt/2. 

In figure [6] we have presented the dependence of the effective azimuthal pressure induced by 
the cylindrical shell on the wedge sides, a 4 p2 Cy i, as a function of r/a for different values of the 
parameter q. 
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Figure 2: The part in the VEV of the energy density, ci 4 (Tq ) cy i, induced by the cylindrical 
boundary as a function on x = (r/a) cos (p and y = (r/a) sin (p for a wedge with (fio = tt/2. 




Figure 3: The part in the VEV of the radial stress, a A (T^} cy \, induced by the cylindrical boundary 
as a function on x = (r/a) cos (j) and y = (r/a) svacj) for a wedge with 0o = tt/2. 
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Figure 4: The part in the VEV of the azimuthal stress, a, (T% ) cy i, induced by the cylindrical 
boundary as a function on x = (r/a) cos (p and y = (r/a) sin (p for a wedge with (fro = tt/2. 




Figure 5: The part in the VEV of the off-diagonal component, a 4 (T 1 2 ) cy i, induced by the cylin- 
drical boundary as a function on x = (r/a) cos <p and y = (r/a) sin (f> for a wedge with <pQ = tt/2. 
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Figure 6: The effective azimuthal pressure induced by the cylindrical shell on the wedge sides, 
a 4 p2cyb a s a function of r/a. The numbers near the curves correspond to the values of the 
parameter q. 

There are several special cases of interest for the geometry of boundaries we have consid- 
ered. The case (f>o = tt corresponds to the semi-circular cylinder. The Casimir energy for the 
corresponding interior region is evaluated in Ref. [34] by using the zeta function technique. The 
case 4>q = 2tt corresponds to the geometry of a cylindrical shell with a coaxial half-plane. And 
finally, the limit (po — * 0, r, a — * oo, assuming that a — r and a<po = b are fixed, corresponds 
to the geometry of two parallel plates separated by a distance b, perpendicularly intersected 
by the third plate. In the latter case it is convenient to introduce rectangular coordinates 
(x , x' 2 , x' 3 ) = (x,y,z) with the relations x = a — r, y = r<p. The components of the tensors 
in these coordinates we will denote by primes. The corresponding vacuum energy-momentum 
tensor is presented in the form 

(0\Tj*\0) = (if)® + (if)® , (56) 

where (T^ )(°) is the vacuum expectation value in the region between two parallel plates located 
at y = and y = a and (Tj?)^ 1 ' i s induced by the intersecting plate at x = 0. The latter is 
related to the quantities investigated above by formulae 

<7f >« = lim <7?> cyl , (#>W = - lim -{T^, (57) 

a 

with lim corresponding to the limit a — * oo, <pQ — » for fixed a — r and acpQ. Taking this limit 
in the term with m = of formula (|38p we replace the modified Bessel functions by the leading 
terms of the corresponding asymptotic formulae for large values of the argument and the integral 
is taken elementary. For the terms with m ^ in formulae (138 j) . (j39|) we note that in the limit 
under consideration one has q = ir/cfro — ► oo, and the order of the Bessel modified functions tends 
to infinity. Introducing a new integration variable x — > qmx, we can replace these functions by 
their uniform asymptotic expansions for large values of the order. After these replacements the 
integration and the further summation over m are done by using the formulae from [46j . 

4 Vacuum densities in the exterior region 

In this section we consider the VEVs for the field square and the energy-momentum tensor in the 
region outside the cylindrical boundary (region II in figure H]). The corresponding eigenfunctions 
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for the vector potential are obtained from formulae ^) by the replacement 

Jqmfrr) -> ffjm (7a, TT) = Jqm{ir)Y q ^(^a) - Yqm(jr)J^(ja), (58) 

where, as before, A = 0, 1 correspond to the waves of the electric and magnetic types, respec- 
tively. Now, the eigenvalues for 7 are continuous and in the normalization condition Q the 
corresponding part on the right is presented by the delta function. As the normalization inte- 
gral diverges for 7' = 7, the main contribution into the integral comes from large values of r and 
we can replace the cylindrical functions with the argument jr by their asymptotic expressions 
for large values of the argument. By this way it can be seen that the normalization coefficient 
in the exterior region is determined by the relation 

= y 6 ™^ [4- )2 (7«) + 4" )2 (7«)] • (59) 

Substituting the eigenfunctions into the corresponding mode-sum formula, for the VEV of the 
field square one finds 



m=0 J ~°° Ji} A=0,1 



7 3 J7 ( ^[*SW.^(70,7r)] (6Q) 



\A 2 + t 2 4t (7«) + *#r (70) 



where the functions g&rx) $^(0)) 5gm(7 a > 7 r ) are defined by relations Q, (fTUjl with the 

function f{x) = gqm{l a i x). To extract from this VEV the part induced by the cylindrical shell, 
we subtract from the right-hand side the corresponding expression for the wedge without the 
cylindrical boundary. The latter is given by formula (|14|) . The corresponding difference can be 
further evaluated by using the identity 



1 2 

-\Y, 



ll? L «"")[*SW.ffS(T')]. (61) 



s=l 



where H q m '(z) are the Hankel functions. In order to transform the integral over 7 with the 
last term on the right of (|6ip . in the complex plane 7 we rotate the integration contour by the 
angle tt/2 for the term with s = 1 and by the angle —tt/2 for the term with s = 2. Due to the 
well-known properties of the Hankel functions the integrals over the corresponding parts of the 
circles of large radius in the upper and lower half-planes vanish. After introducing the modified 
Bessel functions and integrating over k with the help of formula (|26p . we can write the VEVs of 
the field square in the form (|13() . where the part induced by the cylindrical shell is given by the 
formula 

^ 2 )cyi = IE' E r^x 3 ^|^G^)[$W(0),^ m (xr)]. (62) 

m=0 A=0,1 K qm \XO) 

In this formula the functions G^ r,FX ^ [<J?(<^>), /(x)] are defined by formulae (|28p. (|29p . Comparing 
this result with formula (j27[) . we see that the expressions for the shell- induced parts in the 
interior and exterior regions are related by the interchange I qm <^ K qm . The VEV ([62]) diverges 
on the cylindrical shell with the leading term being the same as that for the interior region. At 
large distances from the cylindrical shell we introduce a new integration variable y = xr and 
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expand the integrand over a/r. For q > 1 the main contribution comes from the lowest mode 
m = and up to the leading order we have 

For q < 1 the dominant contribution into the VEVs at large distances is due to the mode m = 1 
with the leading term 



COS ( 2 ^) , ( _ 1) VFll±l_ 

2q + 3 1 ' 2g+l 



(64) 



For the case q = 1 the contributions of the modes m = and m = 1 are of the same order 
and the corresponding leading terms are obtained by summing these contributions. The latter 
are given by the right-hand sides of formulae ([63]) and ([64"]) . As we see, at large distances the 
part induced by the cylindrical shell is suppressed with respect to the part corresponding to the 
wedge without the shell by the factor (a/r) 2/3 with f3 = min(l,g). 

Now we turn to the VEVs of the energy-momentum tensor in the exterior region. Substi- 
tuting the eigenfunctions into the corresponding mode-sum formula, one finds (no summation 
over i) 



^ ^ o J ^ y^r^ jgm ( 7 «) + y s 



/«[*$(0),S$(7o,7r)] 



qm 2 {ia) 



(65) 



(01^10) = rnsm&irn*) <ttjT ^7 E C-^ 1 ^^ - (66) 



m=0 """" " u A=0,1 



Subtracting from these VEVs the corresponding expression for the wedge without the cylindrical 
boundary, analogously to the case of the field square, it can be seen that the VEVs are presented 
in the form (|37p . with the parts induced by the cylindrical shell given by the formulae (no 
summation over i) 

(^)cyl = /sE'E / ^x 3 ^^FW[cI>W(0),K gm (xr)], (67) 
n m=0 A=0,1 K qm (xa) 

(T 2 % yl = f^f E msin ( 2 ^) E(- 1 ) A / dxx^^K^(xr). (68) 
47r ar m=0 a=o,i - 70 Kfri(xa) 

Here the functions -F^- 1 [&{<t>), f{y)] are defined by formulae (|40p . (14ip . By using the inequality 
given in the paragraph after formula (03]), we can show that the vacuum energy density induced 
by the cylindrical shell in the exterior region is positive. 

In the way similar to that for the interior region, for the force acting on the wedge sides 
is presented in the form of the sum (|46p . where the part corresponding to the wedge without 
a cylindrical shell is determined by formula (|47p and for the part due to the presence of the 
cylindrical shell we have 

q °° i f°° I^m (xa) 

P2cyl = -(r|)cylU=0^o = --fE E / dx ^ ^ - F qm[ K qm{xr)]. (69) 

m=0 A=0,1 K qm (xa) 

In this formula, the function F„ [f(y)] is defined by relations (|49p and the corresponding forces 
are always attractive. 
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The leading divergence in the boundary induced part (|67|) on the cylindrical surface is given 
by the same formulae as for the interior region. For large distances from the shell and for q > 1 
the main contribution into the VEVs of the diagonal components comes from the m = 0, A = 1 
term and one has (no summation over i) 



i T i)cy\ 



qci 



15vr 2 r 4 



Co = c 3 = 2, ci = 1, c 2 



-5. 



(70) 



In the case q < 1 the main contribution into the VEVs of the diagonal components at large 
distances from the cylindrical shell comes from the m = 1 mode. The leading terms in the 
corresponding asymptotic expansions are given by the formulae 



(^) cyl «-g 2 (g + lh(g) 



cos(2g0) (a\ 2 i 



7r 2 r 4 



with the notations 



co(q) = c 3 (q) 



2q + 3 



, ci(<?) 



2q 2 + q + 1 
(2q + l)(2q + 3) 



C2(q) 



q+l 
'2q+l 



(71) 



(72) 



In the case q = 1 the asymptotic terms are determined by the sum of the contributions coming 
from the modes m = and m = 1. The latter are given by formulae (|70p . (|7ip . For the off- 
diagonal component, for all values q the main contribution at large distances comes from the 
m = 1 mode with the leading term 



( T 2>cyl 



q 3 (q + 1) sin(2c/(/>) /a\ 2 9 
2q + 1 7r 2 r 3 Vr / 



(73) 



For large values of q, q S> 1, the contribution of the terms with m > is suppressed by the factor 
exp[— 2qmln(r / a)} and the main contribution comes form the m = term with the behavior 
(F 2 ) cy \ oc (/ and (Tf )cyi °c g. In figure[7]we have plotted the dependence of the effective azimuthal 
pressure induced by the cylindrical shell on the wedge sides, a 4 p2cyi> as a function of r/a for 
q = l. 



-0.2 



■0.4 



-0.6 



-0.8 




Figure 7: The effective azimuthal pressure induced by the cylindrical shell on the wedge sides, 
a 4 p2cyb as a function of r/a in the exterior region for q = 1. The curves for the values q = 0.5, 3 
are close to the plotted one. 
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5 Conclusion 



In this paper we have investigated the polarization of the electromagnetic vacuum by a wedge 
with a coaxial cylindrical boundary, assuming that all boundaries are perfectly conducting. Both 
regions inside and outside of the cylindrical shell (regions I and II in figure [1]) are considered. In 
section[2]we have evaluated the VEVs of the field square in the interior region. The corresponding 
mode-sums contain series over the zeros of the Bessel function for TM modes and its derivative 
for TE modes. For the summation of these series we used a variant of the generalized Abel- 
Plana formula. The latter enables us to extract from the VEVs the parts corresponding to the 
geometry of a wedge without a cylindrical shell and to present the parts induced by the shell 
in terms of integrals which are exponentially convergent for points away from the boundaries. 
For the wedge without the cylindrical shell the VEVs of the field square are presented in the 
form (|23p . The first term on the right of this formula corresponds to the VEVs for the geometry 
of a cosmic string with the angle deficit 2tt — <f>Q. The angle-dependent parts in the VEVs of 
the electric and magnetic fields have opposite signs and are cancelled in the evaluation of the 
vacuum energy density. The parts induced by the cylindrical shell are presented in the form 
(I27j) . We have discussed this general formula in various asymptotic regions of the parameters 
including the points near the edges and near the shell. In section [3] we consider the VEV of the 
energy-momentum tensor in the region inside the shell. As for the field square, the application 
of the Abel-Plana formula allows us to present this VEV in the form of the sum of purely wedge 
and shell- induced parts, formula (I37h . For the geometry of a wedge without the cylindrical 
boundary the vacuum energy-momentum tensor does not depend on the angle <j) and is the same 
as in the geometry of the cosmic string and is given by formula (I45p . The corresponding vacuum 
forces acting on the wedge sides are attractive for (f>o < tt and repulsive for <f>Q > n. In particular, 
the equilibrium position corresponding to the geometry of a single plate (4>o = tt) is unstable. 
For the region inside the shell the part in the VEV of the energy-momentum tensor induced 
by the presence of the cylindrical shell is non-diagonal and the corresponding components are 
given by formulae (f38l) . (i39j) . The vacuum energy density induced by the cylindrical shell in the 
interior region is negative. We have investigated the vacuum densities induced by the cylindrical 
shell in various asymptotic regions of the parameters. For points near the cylindrical shell the 
leading terms in the asymptotic expansions over the distance from the shell are given by formulae 
(I50p . These terms are the same as those for a cylindrical shell when the wedge is absent. For 
a wedge with <J)q < tt the part in the vacuum energy-momentum tensor induced by the shell 
is finite on the edge r = 0. For (j)Q > tt the shell-induced parts in the energy density and the 
axial stress remain finite, whereas the radial and azimuthal stresses diverge as r 2 ^ 71 1 '^o -1 ). The 
corresponding off-diagonal component behaves like 7- 27r /</>o-i f or a ji values <f>o- For the points 
near the edges (r = a, 4> = 0, (fro) the leading terms in the corresponding asymptotic expansions 
are the same as for the geometry of a wedge with the opening angle </>o = vr/2. In the limit of 
small opening angles, 4>o <C tt, the shell-induced parts behave like l/0o- I n the same limit the 
parts corresponding to the wedge without the shell behave as 1/<Pq, and for points not too close 
to the shell these parts dominate in the VEV of the energy-momentum tensor. The presence 
of the shell leads to additional forces acting on the wedge sides. The corresponding effective 
azimuthal pressure is given by formula (|48p and these forces are always attractive. 

The VEVs of the field square and the energy-momentum tensor in the region outside the 
cylindrical shell are investigated in section [H As in the case of the interior region, these VEVs 
are presented as sums of the parts corresponding to the wedge without the cylindrical shell and 
the parts induced by the shell. The latter are given by formula (|62p for the field square and by 
formulae (|67p . (|66|) for the components of the energy-momentum tensor. In the exterior region 
the vacuum energy density induced by the cylindrical shell is always positive. Additional forces 
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acting on the wedge sides due to the presence of the shell are given by formula (|69|) . As in the 
case of the interior region these forces are attractive. For large values of the parameter q, the 
contribution into the parts induced by the cylindrical shell coming from the modes with m/0 
is exponentially suppressed, whereas the contribution of the lowest mode m = is proportional 
to q. Though in this limit the vacuum densities are large, due to the factor 1/q in the spatial 
volume element, the corresponding global quantities tend to finite limiting values. 
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